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I. Introduction 

Since 1975 members of the Department of Geography, Columbia University, 
have worked with NASA scientists at the Goddard Institute for Space Studies 
(GISS) on research concerned with satellite remote sensor observations of 
earth resources. This activity is currently supported under NASA Cooperative 
Agreement NCC 5-20. During 1981 research activities have concentrated in four 
areas ; 

1) Simulation of Landsat D Thematic Mapper (TM) observations from aircraft 
multi spectral scanner data and field spectrometer data collected over a 
corn-soybeans agricultural region in Webster County, Iowa during the 1979 
growing season in support of the NASA/AgRISTARS program; 

2) Analysis of the simulations to evaluate the potential utility of the TM 
band 5 (1.55 -1.75 pm) mid-infrared observations in corn-soybeans 
discrimination; 

3) Analysis of current Landsat data to study snow cover in northern New 
England and wetlands in Nebraska and Vermont in cooperation with scientists at 
the U.S. Army Corps of Engineers Cold Regions Research Laboratory; 

4) Application of satellite remote sensor data in selected additional 
agro-environmental research areas pertinent to Columbia staff interests. 

Dr. Samuel N. Goward, v'esearch associate, and Professor Robert Lewis, 
chairman of the Geography Department, are the principal investigators. Dr. 
Goward serves as the research director. Professor Leonard Zobler supports the 
activity as faculty advisor. Columbia University students work on projects 
under investigation as research assistants and aides. During 1981 five 
graduate students were funded under the cooperative agreement. Several other 
students, funded from other sources, also supported the research activities. 
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The research activities are carried out at the Goddard Institute for Space 
Studies. The Institute provides work space, data, computer facilities, and 
software to support the research. The Columbia staff interact daily with the 
NASA scientists to carry out the Earth Resources Group research activities. 
Principal GISS contacts are Dr. Stephen Ungar and Dr. Richard Kiang. Senior 
Columbia staff collaborate with the NASA scientists to plan and direct the 
research activities. The research assistants and aides participate in the 
research and provide technical support in preparation of presentations and 
reports. 

II. Research Activities 

A. Simulation of Landsat D Thematic Mapper Observations 

In the third quarter of 1982, NASA plans to launch Landsat D. A new 
observing system, the Thematic Mapper, will bs flovrn on this mission. Thematic 
Mapper represents a significant technical advance compared to the current 
Landsat observing system (J.L. Engel, 1980). A comparison of relevant TM and 
MSS parameters is given in Table 1. Considerable prelaunch research activity 
is underway to evaluate the potential information gain TM data will bring to 
earth resources observations (ORI, 1981). Columbia/GISS Earth Resources 
scientists are contributing to these activities through analysis of field 
measurements data collected at AgRISTARS "supersites". 

The AgRISTARS field measurements program collects Landsat and ground 
observations at 300 selected agricultural sites across the United States 
(AgRISTARS, 1981). In addition, at two "supersites" (Webster County, Iowa for 
corn and soybeans and Cass County, North Dakota for small grains) more 
intensive observations are collected by the helicopter-based field 
spectrometer system (FSS), the NSOOl aircraft-based multispectral scanner 
system and by ground observers. (For a discussion of the analogous LACIE field 
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TABLE 1 


Comparison of Landsat MSS and Thematic Mapper 
Sensor Specifications 


LANDSAT THEMATIC MAPPER 



Band 

H 

Spectral 

Band # 

Spectral 




Interval (pm) 


Interval (ym) 


1 


0.50-0.60 

1 

0.45-0.52 


2 


0.60-0.70 

2 

0.52-0.62 

Spectral 

3 


0.70-0.80 

3 

0.63-0.69 

Resolution 

4 


0.80-1.10 

4 

0.76-0.90 





5 

1.55-1.75 





6 

2.08-2.35 





7 

10.40-12.50 

Spatial 


80 

meters 

30 

meters 

Resol uti on 




120 

meters Band 7 

Radi ometri c 



26 


2® 

Resol uti on 

(64 

gray levels) 

(256 gray levels) 
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measurements program see Bauer, et a1. 1978.) The supersites consist of 5 x 6 

nautical mile regions in which 80 fields are periodically observed by field 
enumerators (NASA-USDA/ESS, 1979). These observations are timed to coincide 
with Landsat overpasses and are conducted at a 9 to 18 day interval. The 
helicopter and aircraft observations are also conducted at the same time 
although the aircraft is flown less frequently (once in 1979, 4 times in 1980 
and 1981). The aircraft and Landsat observations provide multi spectral 
scanner data for the entire site. The helicopter FSS observations are 
non-imaging high spectral resolution radiance measurements from a flight path 
across each of the 80 periodically observed fields. 

The Conceptual Framework of the GISS Simulations 

During 1981, the Columbia/GISS Earth Resources Group, as part of their 
AgRISTARS research effort, have developed techniques to simulate TM 
observations from the field measurements data. The satellite, aircraft, and 
helicopter observations serve complimentary roles in the TM simulation. As 
shown in Figure 1, the aircraft data serve as a means to simulate both TM 
and Landsat observations. The Landsat simulation may be compared to the real 
Landsat data to confirm the simulation technique. The FSS data serve as a 
calibration standard to verify radiometric adjustment procedures, and as the 
basis for extending the aircraft observations to satellite observing 
conditions through simulation of atmospheric effects and consideration of the 
satellite TM engineering parameters. In addition, the FSS observations permit 
temporal evaluation of crop radiance behavior in TM bands at those times when 
the aircraft is not flown. 

The GISS-developed techniques include radiometric, geometric, and 
resolution processing of the aircraft scanner data and radiometric and 
statistical processing of the FSS data. These techniques were tested using 
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FIGURE 1 TM SIMULATION APPROACH ' j 

! 

From Presentation material used by S. Ungar; Renewable Resources Thematic I 

Mapper Simulator Workshop (ORI, 1981). I 
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1979 field measurements data from the Webster County, Iowa supersite (Figure 
2). One NSOOl flight, collected on 8/30/79 and four FSS observations 
collected on 6/11, 6/29, 7/16, and 8/30/79, were processed to simulate TM 
observations. 

TM Simulations From NSOOl Data 

The NSOOl (NASA scanner 1) instrument was constructed by NASA engineers at 
the Johnson Space Center to emulate, as much as possible, the satellite TM 
instrument (Richard, Merkel and Meeks, 1978). The NSOOl instrument parameters 
are given in Table 2. The NSOOl is typically flown on the NASA C-130 aircraft 
at altitudes between 5 and 8 kilometers. With a 2.5 milliradian detector 
IFOV, the scanner observes, at nadir, a 19 meter square ground resolution 
element from 7,620 meters (25,000 feet) altitude. At this altitude, the 
scanner normally operates at 15 revolutions per second which, for the nomimal 
aircraft ground speed of 135 meters/sec, produces a 50% overlap between 
adjacent scan lines. As the instrument was constructed prior to final 
selection of TM bands, eight rather than seven bands are included in the 
detector package (NSOOl band 5, 1.00 - 1.30 pm is the non-TM band). 

Although the NSOOl scanner was constructed to emulate the satellite TM 
scanner, no aircraft scanner system is fully capable of simulating satellite 
observing conditions. Not only does the aircraft fly at much lower altitudes, 
thus not observing the full effects of atmospheric attenuation, but also, 
because of the altitude difference, the scanner must observe the surface at 
more extreme scan angles to view a ground area comparable to a small subarea 
in one a.'fellite scene (OR I, 1981). 

The NSOOl maximum scan angle from nadir is 50 degrees. This permits the 
instrument to view an across track ground swath of 18 kilometers (± 9 
kilometers from nadir) at an altitude of 7,620 meters. This variation in look 
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TABLE 2 

NSOOl TECHNICAL SPECIFICATIONS 


srrcrRA:L bands 


Ban'l Detector Spectral bandwidth, ■■ m NET. » 


1 

S; 

0,45 

- 

0.52 

0,5 

2 

S. 

a 

0,52 

- 

0,60 

.5 

3 


0.63 

- 

0.69 

.5 

4 

*• 


0,76 

- 

0.90 

.5 

5 

r-e 

1.00 

- 

1,30 

1.0 

f 

ve 

l.bS 

- 

1 , 7.5 

1.0 

7 

In^s 

2.06 

• 

2.21 

2.0 

B 

HgCdTe 

10,4 

- 

12,50 

NE,'.T 


DJS’.CN DATA 


Instantaneous field of view (IFOV) 
AcrosE-tracli field of view 
Nor'inal aperture diameter 
Effective aperture area 
( numtier 

Frir.ary focal length 
Inflight calibration 

Short wavelength array temperature 
V ‘H ranee 
Scan rate 

Scan speed stability 

Data quant 1 cation 

N-rber of video sa~ples. scan line 

p:ll compensation 

Data format 

Soar, mirror 

Tape packing density 

Recording code 


2.5 * 2. S mill iradians 
100 ' 

10,16 cm 
72,4 cm^ 

1.15 
16.6 cm 

Integrating aphere and two 
controllable blackbodiea 

19B'K 

Variable 0,025 to 0.25 

Variable 10 to 100 acans/sec 

One'- third of the IPOV, acan line to 
scan line 

6 bits (256 discrete levels) 

700 

;15' 

Com.patible with M^S 

45* rotating mirror 

10 000 biti-per inch (BPI) constant 

Bi-«-L 
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angle introduces significant across track variations in ground resolution and 
the radiometric conditions observed by the sensor that will not be present in 
the TM observations. 

An additional factor, which should not be present under satellite 
conditions, is complex distortion of image geometry due to short term changes 
in aircraft motion. At 8,000 meters altitude, the C-130 aircraft is subject to 
uncontrolled motion due to buffeting in the turbulent lower atmosphere. Short 
term changes in the velocity, altitude, pitch, yaw and roll of the aircraft 
can significantly alter the geometry of the image data in comparison to 
satellite observations (and other information available such as photography 
and maps). The NSOOl is roll -compensated (Richard, Merkel, and Meeks, 1978) 
Hov^ever^, other motions of the aircraft during data acquisition significantly 
affect the image geometry. Figure 3 presents a grayscale image of raw NSOOl 
Band 6 (TM5) data from 8/30/79 v/hich displays the resolution, radiometric and 
geometric prperties of a typical aircraft scanner image. 

Aircraft Scanner Data Processing Techniques Developed 

Although processing techniques developed by the Columbia/GISS staff were 
specifically directed to processing the NSOOl data, these same techniques are 
applicable to aircraft scanner data from other systems such as the NASA/ERL 
TM simulator and the NASA/Ames Daedal as system. The approach is to develop 
techniques, where possible, that systematically account data variations 
introduced by the scanner system. For example, ground resolution variations 
are, in general, an explicit function of scan angle. However, certain 
characteristics of the data must be treated empirically. Variations in 
atmospheric optical depth, scattering, and surface bidirectional reflectance, 
although they appear to be related to scan Took angle, are not sufficiently 
known to be subject to systematic correction. In this case, empirical 
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techn1que$ were developed to compensate or adjust the data in a reasonable 
fashion. Correction of the effects of platform velocity, altitude and attitude 
variations on Image geometry presents a particularly difficult problem. 
Although systematic techniques are possible, practical Implementation requires 
a high quality source of aircraft navigational data during scanner 
observations. Under relatively stable flight conditions, empirical 
“rubbersheetlng" techniques appear to be adequate. However, under unstable 
flight conditions a more systematic approach Is needed. Expermental activities 
In both areas were explored by the Columbla/GISS group In 1981. 

Ground Resolution 

Whereas the NSOOl ground resolution Is 19 meters at nadir from a 7,620 
meter altitude, at the 50 degree maximum look angle, the ground resolution 
element Is approximately 30 meters In the along track direction by 46 meters 
In the across track direction. Since the NSOOl Instrument Is roll -compensated 
and pixels are contiguous in the scan direction, calculation of the scan angle 
position of each pixel in a scan Hne Is straight forward. The central pixel 
In each scan line represents the nadir pixel, each adjacent pixel Is 2.5 
mil 11 radians of scan angle away from nadir. The ground resolution of any 
pixel in the scan may be calculated by (Saker and Mikhail, 1975): 
flalong track “ (h/cos A)a (1.1) 
llacross track " (h/cos^ A)a (1.2) 

where R - linear resolution 
h s altitude 
A = look angle 

a = angular resolution of the detector. 

To aid 1n later geometric processing, the data are resampled to nominal 10 
meter pixels. In the scan direction resampling 1s conducted using equation 
1.2. Since adjacent scans are 50% overlapped each scan is assumed to 
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represent a 10 meter ditta sample. After geometric correction the 10 meter 
pixels are aggregated to 30 meter pixels to simulate the nominal satellite TM 
ground resolution. Figure 4 Is a grayscale of the 8/30/79 NSOOl band 6 (TM 
band 5) data which has been resampled to 10 meters and radiometrically 
adjusted. 

Radiometric Adjustment 

Because the NSOOl makes observations to 50 degrees either side of nadir 
the radiance measured varies not 'only as a function of the ground cover 
present but also due to the bidirectional radiance properties of the ground 
cover, variations in atmospheric path length and anisotrophic scattering of 
the radiant flux by atmospheric aerosols. The relation between incident solar 
radiation and the NSOOl maximum look angles for the 8/30/79 observation is 
plotted in Figure 5. Note that the solar flux is nearly perpendicular to the 
eastern look angle and nearly parallel to the western look angle. Figure 6 is 
a plot of the mean and standard deviations for each column of the western 
portion of the 8/30/79 NSOOl band 1 (0.45 - 0.52 urn) data. The general 
pattern displayed is increasing radiance with increased westerly look angle; 
less significant variations occur in the easterly direction. The reduced canopy 
shadowing and increased canopy cover along with increased observed aerosol 
scattering due to both increased path length and increased back scatter 
parallel to incident solar illumination, cause the image to brighten in the 
westerly direction. 

The Column-Averaging Approach 

The landscape and atmospheric conditions that contribute to this effect 
are not sufficiently known to explicitly compensate for their variation as a 
function of look-angle. A simple empirical approach to radiometric 
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adjustment, which will insure a uniform tonal distribution across the image, 
is to derive a least-squares fit to the column-average distribution. As shown 
in Figure 6, this least-squares line can be used to calculate a 
column-by-column adjustment of the data that produces the uniform tonal 
appearance of the image data displayed in Figure 4. 

The assumptions that underlie this column-average radiometric adjustment 
approach are that landscape materials of differing reflectance properties are 
randomly distributed in the scan direction and that the bidirectional 
reflectance of the landscape is independent of landscape type (e.g., corn and 
soybean canopies). From visual inspection of the 8/30/79 NSOOl data the 
former assumption is reasonably satisfied for this data set. The latter 
assumption was statistically examined for the 8/30/779 data. 

Following the classification procedures discussed later in this report, 
the radiometrical ly adjusted data were classified, on 0^ commission error 
basis, to identify corn and soybeans pixels. The classification map was used 
as a mask on the raw data to select corn and soybeans pixels. Data 
column-averages and least-squares fits were then computed independently for 
corn and soybeans. Any difference in the corn and soybeans column average 
patterns would suggest differential corn/soybeans bidirectional reflectance as 
a function of look angle. Visual inspection of the comparable corn/soybeans 
plots revealed little difference in the two curves. Statistical evaluation of 
the least-squares fits indicated that there is a slight but statistically 
significant difference. Thus, although differences in corn and soybean 
canopies during this 8/30/79 observation appear to differentially effect the 
column-average pattern of the data, the difference is small and can be, for 
this data set, ignored. 
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Comparison of Ra'4iometrica1ly Adjusted NSOOl Data and FSS 
Observations 

The FSS observations may be used as an independent check of the 
radiometric adjustment procedure. The FSS instrument acquires nadir 
observations of the same ground area the NSOOl scans. FSS data are processed 
to NSOOl sprectral band intervals (see "TM Simulations from FSS Data"). By 
comparing the two data sets before and after NSOOl data radiometric adjustment 
the effectiveness of the adjustment may be evaluated. 

The locations of the helicopter flight lines were derived by analysis of the 
FSS borcslght color photography and C'lR photography acquired during the NSOOl 
overflight. The FSS instrument, when flown at 60 meters, observes a nominal 
24 meter spot size (Bauer et al., 1978). However, one rotation of the filter 
wheel takes about one second, during which time the helicopter moves 
approximately 27 meters with respect to the ground. To select equivalent 
NSOOl observations every other NSOOl pixel along the helicopter flightline was 
extracted from the NSOOl data. The mean value per field of the FSS and NSOOl 
data were correlated to study the impact of the radiometric adjustment 
procedure. 

Table 3 presents the results from the correlation analysis. Correlations 
between the two data sets prior to adjustment are consistently lower than 
following adjustments. This indicates that radiometric adjustment makes the 
NSOOl data appear more like nadir observations. The relation, after adjustment 
for all crops shows a strong relation (R^ > 0.9) between the two data sets. 
Correlation between corn and soybean means taken separ-ately shows somewhat 
lower final R^ values. This is most likely the result of low within-crop 
variability in the reflectance of corn and soybean canopies, close to the 
sensitivity limits of the two instruments. However, the differential 
bidirectional reflectance of the corn and soybeans canopies may also affect 
this relation. Scatterplots for NSOOl bands 1 and 4 of the before and after 
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adjustment relations between FSS and NSOOl data are given in Figure 7. These 
plots confirm the tabular results from the correlation analysis. 

Geometric Correction 

Aircraft altitude, velocity, and attitude variations during scanner 
observations introduce complex geometric distortions in the scanner imagery. 
These distortions will generally not be present in the satillite data and make 
comparison of the aircraft data with maps and photos of the area, as well as 
construction of registered temporal data sets, difficult. Two approaches, 
one empirical and the other systematic, to geometric correction of the 
aircraft scanner data have been investigated by the Columbia/GISS staff during 
1901. 

The empirical geometric correction approach uses map and image control 
points to statistically de>'ive a "rubbersheet" transform of the image data 
that matches, as much as possible, the geometry of the map control point 
source (Baker and Mikhail, 1975). The systematic approach use aircraft 
altitude, velocity, and attitude information, recorded during the observation, 
to compensate for the effects of these variations on the image geometry 
(Spencer, Wolf, and Schall, 1974). Although the latter approach is more 
desirable, since it represents a "true" correction of these distortions rather 
than a statistical estimate, it is highly sensitive to the precision and 
accuracy of the recorded aircraft altitude, velocity, and attitude 
information. The experiments completed at GISS during 1981 suggest that for 
the current NSOOl observations the empirical approach is more sucessful than 
the systematic approach. 
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FIGURE 7 

COMPARISON OF FSS AND NSOOl MEANS BEFORE AND AFTER 
NSOOl RADIOMETRIC ADJUSTMENT 
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Empirical Geometric Correlation 

The empirical "rubbersheet" approach to geometric correction requires a 
dense network of control points that can be easily located in the image data 
and a ground reference source. For the United States, USGS 7.5 minute 
quadrangles provide an accurate source of ground information. In the Webster 
County, Iowa, region one feature, which is common to the maps and and the image 
data, is the road network. The road system follows the Township and Range 
survey system at one mile intervals. The intersections between the N-S/E-W roads 
were selected as ground control points. The road intersections are uniformly 
distributed in the data and easily located in the maps and in the image data. 

To encompass Webster County supersite 59 road intersections were selected as 
ground control points (figure 8). 

The map coordinate locations of the road intersections are derived from 
the USGS maps by locating the center of each intersection on the map and 
digitizing the locations via the flatbed digitizer at the Columbia University 
Lamont“Dorherty Laboratories. The digitizer table precision is 1/1000 of an 
inch. For the USGS 1/24,000 scale map this results in a ground coordinate 
location precision uf better than one meter. 

The image coordinates ^are located on line printer grayscales of the ten 
meter data. Preliminary experiments were conducted using grayscales of NSOOl 
bands 3,4, and 6 to determine which band or bands were best suited for 
identification of the road intersections. Band 3 was found to be most 
effective, followed by band 6 and band 4. Two analysts, each working with a 
different bands, independently identified the line and column locations of the 
59 road intersections in 10 meter data. The road intersection is defined as 
the point at which two lines, one located at the center of the N-S road and 

one located at the center of the E-W road, intersect in the grayscale (Figure 
9). For the 8/30/79 data the two analysts generally selected the same ten 
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ineter pixel. In those cases where they disagreed the difference was one 
coordinate location in one direction only. 

The rubbersheet transform is derived by calculating a least-squares 
polynomial relation between the image and map coordinates. The general form of 
the relation is; 

1' “ ao ai x< (E.l) 

y' ■ bo + j b, (2.2) 

1*1 

where x',y' represent the map coordinates and x,y represent the image 
coordinates. This relation is used to reassign data from the image coordinate 
system to the map coordinate system. The nearest neighbor criterion is used 
to specify which image data are assigned to each map coordinate location. 

Transformation experiments conducted with the 8/30/79 NSOOl data led to 
the conclusion that a third order polynomial fit produces visually the most 
acceptable solution. Although higher order polynomials result in lower root 
mean square (RMS) error relations between the image and map control points 
than the third order, they also introduce geometric distortions in the 
resultant image that are not present in the raw data. The RMS error to a 
third order fit on the 59 control points in the 3/30/79 data is approximately 
60 meters (two TM pixels). 

Local Rubbersheeting 

Analysis of the third-order transformed 8/30/79 data Indicated that 
the greatest residual error was a high frequency variation of the distance 
between control points in the flightline direction (Figure 10); possibly due 
to aircraft pitch motion. Higher order polynomials had been judged 
ineffective. An alternate approach in this case is to rubbersheet local 
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strips of the data in which this residual effect could be treated by lower 
order polynomials (Figure 11). To insure continuity of the image strips 
following transformation, the transformation of each strip is constrained to 
fit selected control points in adjacent strips. This insures that the 
transformed image strips will match at these points and, if the local 
transformation is reasonable, that the strips will properly mosaic. Local 
rubbersheeting of the 8/30/79 data produced a visually acceptable image with 
an average RMS error for the 59 control points of less than 10 meters (Figure 
12 ). 

Systematic Gtometric Correction 

Experiments with the 8/30/79 NSOOl data demonstrated that the empirical 
rubbersheet approach to geometric correction is not well suited to aircraft 
scanner images v^here short term variations in platform motion have introduced 
local geometric distortions. Preliminary screening of the 1980 AgRISTARS NSOOl 
observations indicated that this high frequency local geometric distortion is 
typically present in the NSOOl data. The 9/10/80 flight over Webster County, 

Iowa presents an extreme case of this problem (Figure 13). There is a great 
need to establish a more systematic approach to geometric correction o^ 
aircraft scanner data which draws upon recorded aircraft motion information to 
correct for motion-introduced distortions. 

NERDAS 

The navigational parameters of the NASA C-130 aircraft are recorded 
during observation missions by the NERDAS (NASA Earth Resources Data 

> ■ 

Annotation System). Signals from the Litton 51 navigational computer [ 

(LTN-51), radar altimeter, and attitude sensors (pitch and roll) as well as | 

I 

clock time are recorded by the NERDAS (NASA, 1980). The relevant NERDAS 
recorded parameters and their precision is given in Table 4. 
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TABLE 4 


NERDAS-RECORDED NAVIGATIONAL PARAMETERS NEEDED 
FOR SYSTEMATIC GEOMETRIC CORRECTION 


PARAMETER 

Time 

Altitude 

Heading 

Drift 

Roll 

Pitch 


NERDAS PRECISION 


0,1 second 
10 feet 
0.1" 

0 . 1 " 

0 . 1 " 


0 . 1 ° 


Ground Speed 


1 knot 



In the ideal case the NERDAS data should be updated for each observation 
(pixel) of each scan. At a s>can rate of 15 revolutions per second (rps) an 
update rate of 0.2 milliseconds would be required. Under the assumption that 
the aircraft is stable during one scan^ an update rate of 67 milliseconds is 
needed at 15 rps to account for motions between scans. The NERDAS records are 
updated every 40 milliseconds, however individual parameters have differing 
update rates (Spencor, Wolfe, and Schall, 1974). In particular, information 
from the LTN-51, including ground speed, heading, and drift, are updated every 
1.8 seconds. This is convoluted with the 1 second LTN buffer update rate 
which creates at least a 0.5 second uncertainty in what time the measurements 
were taken. At the 15 rps scan rate 27 .t 7.5 scans will occur between updates 
of these critical motion parameters. 

The equations employed to systematically correct the geometry of the 
NSOOl data are given in Table 5. Note that since the NSOOl instrument is 
roll -compensated this parameter is not included in the equations. Clock time 
recorded in the NERDAS and NSOOl data are used to register the data sets. The 
LTN-51 derived parameters are cubicly interpolated and linearly smoothed to 
reduce the discontinuous nature of the observations. Bad NERDAS records, as 
noted in the OSC NERDAS data quality report, are replaced with interpolated 
values. 

Preliminary experiments employing the NERDAS data to systematically 
correct the NSOOl image geometry have been dissappointing. The resultant 
image for the 8/30/79 data is shown in Figure 14. Analysis of this image 
compared to ground conditions indicates that the technique and/or NERDAS data 
do not properly account for aircraft variations in drift and/or ground speed. 
Although the LTN-51 data appear suspect, numerous factors, such as the 


TABLE 5 


SYSTEMATIC GEOMETRIC CORRECTION 
USING NERDAS DATA 


original page is 

OF POOR QUALITY 


FOR ROLL COMPENSATED SCAN MIRROR VELOCITY 


X=Xo*h tan0 sin<f'h(1an (j//cos^)cos <j) 
Y=Yo-h tan0 cosi/>-h(tQn if//cos0)sin (j> 


Xo{l*At)*Xo(th/ V sin{(/»*8)dl 

Yo(t*A1)=Yo(t) */’v'cos((/)«B)dt 

where: Xo,Yo= nadir coordinates 
At ' scan period 

0= pitch 
(^'heading 
8= drift 


t|^*scon angle 
h = altitude 

v = ground speed 
f = frequency of 
observation 


UPDATING Xo,Yo WITH NAVIGATION DATA (NERDAS) 

(1) Find nearest observation times from NERDAS 

where; i = lNT {(t-O/f+l) 

(2) Obtain values for each scan by cubic Interpolation 

z(t)= i a„r 

n>0 

where; 2 represents [0,(f),8,;i',h, or v] 
using NERDAS values for Z(tJ i'l = k*i+2 
solve for 0 , n=0,... 4 

(3) Calculate nadir displacement between scans assuming 

rzdl=i 

( n*l I. ' 

where: Z represents v sin((/)+8)or -v cos ((^*8) 
end a„ is determined by technique used in step (2) 
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alignment between the NSOOl and the various motion detectors or Innaccurate 
treatment of velocity effects are other possible sources of error. 

Further analysis of these problems Is required, the systematic approach 
to geometric correction 1s the best general solution to rectification of 
aircraft scanner imagery. The Columbla/GISS staff will explore a number of 
alternatives during 1982. One approach to resolving the current discrepancies 
In the systematic technique will be to examine how well that data can be 
brought into accord with the control points by varying the magnitude of 
individual aircraft motion parameters. 

TM Simulation 

Once the ground resolution, radlometry, and geometry of the NSOOl data 
are processed a simulation of TM observations may be constructed. NSOOl band 
5 Is deleted from the data and the remaining bands, with the exception of 
NSOOl band 8, are aggregated to the TM 30 meter ground resolution. NSOOl band 
B (10.4 - 12.5 pm) is aggregated to 120 meters to simulate the TM thermal 
Infrared band ground resolution. On the basis of the TM orbital parameters the 
data array is rotated and an array of 366 by 313 pixels , representing the 5 
by 6 nautical mile area of the AgRISTARS supersite, Is extracted from the data 
(Figure 15). The data are then convoluted with appropriate header information, 
in JSC Universal format, and forwarded to AgRISTARS scientists in Houston as 
well as subjected to analysis at 6ISS. 

Further steps in processing are possible but have not been carried out in 
1981. Examples include calibration of the NSOOl data to bidirectional 
reflectance factor by regression to the F'^S data; conversion to nominal 
radiance and inclusion of expected atmospheric attenuation at satillite 
altitude; and inclusion of the TM engineering instrument parameters. These 
additional data processing steps will increase the realism of the simulation 
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of the expected satellite TM observations* Further progress in this area Is 
expected In 1982. 

Simulation of TM Observations from FSS Data 

The Field Spectrometer System (FSS), a high spectral resolution (0.02 uni, 
0,4-1. 1 pm; 0,05 ym for 1.1-2. 4 yinj 0.5 pm, 8.0 -14.0 pni) nonimaging 
instrument, periodically observes the AgRISTARS supersites from a helicopter 
platform at a nominal altitude of 60 meters. Missions are typically flown on 
a 9 to 18 day cycUi during the growing season, in phase, as much as possible, 
with Landsat observations of these locations. A set of 10 flightlines are 
flown during each observation period. The flightlines are planned to pass 
over 80 fields that are also periodically observed by USDA field enumerators. 
Instrument calibration is achieved by observing a panel of known reflectance 
before and after each pair of adjacent flightlines are completed. A 
boresighted 70 mm camera is operated in conjunction with the FSS instrument to 
provide a visual assessment of the ground conditions being observed by the 
spectrometer. 

The FSS data are preprocessed at the Johnson Space Center to selectively 
extract only observations of the 80 periodically observed fields and the 
calibration panel. These data are forwarded to LARS-Purdue University, where 
the data are converted to bidirectional reflectance factor by calibration 
to the panel observations. In addition, the field enumerator observations are 
appended to the appropriate FSS observations. The FSS instrument generally 
acquires multiple observations of each field (the number depends on the field 
dimension in the flight direction). Both single-scan and field average data 
are produced by the LARS staff, LARS processing of the field average data 
includes a statistical method to delete "unlike" observations prior to 
computation of the field means (Biehl, 1979). 
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Tht Columb1a/6ISS staff acquire the processed FSS data from LARS and carry 
out further data processing to simulate TM, Landsat and other sensor systems 
spectral configuration. The nominal spectral configuration of the FSS 
instrument is 0.02 um between 0.4 -1.1 0*05 urn between 1.1 -2.4 pm. 

Studies by 6ISS and LARS scientists suggested the actual configuration 
deviates from these specifications. A bench test of the FSS filter wheel 
revealed that the spectral transmission of the filter wheel is a nonlinear 
function of the filter wheel position (Barrnett, 1980). Based on JSC 
recommendations, the recalibrated FSS spectral intervals are used at GISS for 
data processing. TM, Landsat MSS and other spectral reflectances are computed 
from the FSS data by integration of the appropriate FSS observations weighted 
by the proportion of solar illumination in each FSS spectral interval. This 
computation produced the "effective" spectral bidirectional reflectance factor 
that the TM, Landsat MSS or other instrument would observe. Figure 16 
provides an example plot of FSS data and integrated TM band reflectance from 
Webster County, Iowa AgRISTARS supersite. Additional simulation steps that 
may be included to produce more realistic results include; conversion to 
nominal observed radiance, simulation of atmospheric effects observed from 
satellite altitude, inclusion of appropriate sensor filter functions, and 
conversion of radiance measurements to count data based on sensor engineering 
data. 

B. Analysis of TM Simulations 

One of the major reasons to simulate TM satellite observations is to 
provide a data base with which to evaluate the potential information content 
of TM observations for earth resources research. One pressing issue is whether 
the TM band 5 (1.55 -1.75 ym), mid-infrared (or shortwave infrared) radiance 
data will provide information about earth resources that current Landsat MSS 
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bands do not. This band was selected for observation of vegetated landscapes 
on the basis that leaf reflectance variations In these wavelengths are 
principally a function of leaf water conte.it, versus pigments In the visible 
and leaf mesophyll structure In the near-infrared. However, only limited 
research has been conducted to access the comparative utility of the 1.55 - 
1.75 liin band for vegetation analysis. 

In 1981 Columbia/GISS staff conducted a preliminary study of the mId-IR 
(1.55 -1.75 urn) reflectance behavior of corn and soybeans from the data 
collected In Webster County, Iowa. The Investigation Included histogram 
analysis and parallelapiped classification of the processed NSOOl data; a 
comparative analysis of the NSOOl and FSS data for 8/30/79; and temporal 
analysis using the FSS data. The principal conclusion of this work Is that 
the 1.55 -1,75 |im TM band should Improve corn-soybean discrimination In 
comparison to observations restricted to visible and near-infrared portion of 
the spectrum. In particular, for these observations of Webster County, Iowa, 
the 1.55 -1.75 urn measurements provide corn-soybean separability earlier in 
the season than visible and near-infrared observations. 

Experiment Design 

NSOOl data, FSS data and ground periodic observations collected during 
1979 for the Webster County, Iowa AgRISTARS supersite were used for analysis. 
Eighty fields within the site, selected by the field enumerator prior to the 
growing season, are observed every 9 to 13 days; as much as possible, in 
conjunction with other observations by the FSS, NSOOl, and Landsat sensors. 

To produce comparative results, the analysis was constrained to those fields 
which were observed by both the NSOOl and FSS instruments and which could be 
unambigiously located in the NSOOl data. Four fields (29, 30, 37 and 79) were 
not observed by the FSS (apparently due to chickens that were disturbed by the 
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helicopter) and three fields (8,17» and 22) could not be accurately located in 
the NSOOl data due to ambiguities in the available ancillary information. Of 
the 73 remaining fields, 36 were corn, 33 were soybeans and 4 were oats. 
Analysis of the 8/30/79 NSOOl Data 

The NSOOl data were radiometrically adjusted and resampled to 30 meter 
ground resolution, as previously described. The 73 fields were located in 
printer grayscale maps of the data by consultation with AgRISTARS field 
boundary maps and aerial photography taken during the NSOOl overflight and 
earlier in the 1979 growing season. The analysis was constrained to "field- 
center" pixels to study the typical reflectance behavior of corn versus 
soybeans canopies; indepedent of mixed pixels at field borders. Field-center 
pixels were defined by the largest single rectangle which could be located 
within the field without including farmsteads, drainage ditches or other 
non-crop features. 

Classification Approach 

Visual inspection of the NSOOl data on the GISS image display system 
suggested that corn and soybeans could be more easily discriminated in TM band 
5 than in TM bands 2, 3 or 4. To verify this observation, the data were 
subjected to numerical classification using the GISS-ISURSL parallelapiped 
classifier (Hyde, Goward and Mausel, 1977) . 10 corn fields, 10 soybeans 

fields and 1 oats field were randomly selected, on the basis of field number, 
to serve as test fields. The remaining fields served as training fields. 
Several band combinations, specifically bands 2 and 3; bands 2, 3, and 4; 
bands 2, 3 and 6; and bands 2,3,4, and 6, were used to classify the data. 

Preliminary class signatures were extracted from the test field data by 
computing band means and standard deviations for each crop. The parallelapiped 
class signature was then defined, for each band, as the range of data within 
two standard deviations from the mean. Class order is a significant component 
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of the GISS-ISURSL parallelaplped classifier because where class signatures 
overlap a pixel will be assigned to the first class in which it is a potential 
member. The strategy employed is to order class signatures within the 
classifier with respect to signature range magnitude. The class with the 
smallest signature range is entered first. In this case, the corn signature 
was the least variable and the oats signature the most variable. The results 
of training field classification for the preliminary signatures are given in 
Table 6. 

To evaluate the performance of the paraTlelapiped classifier a proportion 
estimate is calculated. This proportion compares the actual number of pixels 
in the class training fields with the numbers of pixels classified as being 
members of that class. Where the proportion is 100, the classified and actual 
are the same and the classification is considered accurate. Proportions 
larger than 100 indicate an overclassification of the category and smaller 
indicate underclassification. Because the parallelapiped classifier is order 
dependent the overclassification of the corn for the bands 2 and 3 combination 
indiacted that the corn signature was too broad and could be impoved by 
reducing the signature range. This signature refinement is conducted 
interactively by reducing the upper and lower bounds of each band until the 
best proportion estimate is achived. Table 7 presents the results of this 
analysis. Note that a significant improvement Is achieved in corn/soybeans 
classification for the bands 2,3 combination. A minor improvement is achieved 
for the 2,3,4 combination and no change was necessary in the remaining band 
combinations. 

Classification Results 

The test field classification results for the 2,3,5 and 2, 3,4,5 TM band 
combinations are of particular interest. In those cases where the TM 1.55 - 
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TABLE 6 

PARALLELAPIPED CLASSIFICATION OF 8/30/79 NSOOl 
PRELIMINARY SIGNATURE 
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TABLE 7 

PARALLELAPIPED CLASSIFICATION OF 8/30/79 NSOOl DATA 
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1.75 Min band is included classification accuracy is improved by approximately 
This suggests that the mid-infrared band will improve our ability to 
discriminate corn and soybeans during this period of the growing season. 
Evaluation of the training fieM performance statistics confirms this 
conclusion (Table 8). Classification accuracy is approximately 6i5 higher with 
TM band 5 included. 

Discussion 

Laboratory investigations of corn and soybeans leaf spectra (LARS, 1968) 
in the 1.55 -1.75 pm spectral range have observed the same differential 
reflectance behavior that is observed in this classification study of NSOOl 
data. In essence, the reflectance behavior of corn and soybeans leaves 
differs most sharply in the 1.55 -1.75 uni portion of the electromagnetic 
spectrum. Little is known about how this leaf reflectance behavior translates 
to canopy reflectance behavior. Examination of histograms of the corn and 
soybeans data extracted from the NSOOl data (Figure 17) suggests that although 
on average the reflectance difference between corn and soybeans about the same 
in the near-infrared and mid-infrared that the variability of the reflectance 
in the mid-infrared is less than in the near-infrared. 

Comparison of 8/30/79 FSS and NSOOl Data 

As previously discussed, the radiometrically adjsuted NSOOl data and the 
FSS data, integrate to TM spectral bands, are highly correlated when field 
means are compared. Further evidence of their equivalence can be noted by 
comparing histograms of individual observations from each data set. The 
regression relation between the NSOOl and FSS field means was used to convert 
the FSS bidirectional reflectance factor data to 8-bit count NSOOl-like 
measurements. Using the FSS single-scan data, histograms of corn and soybeans 
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FIGURE 17 

8/30/79 NSOOl CORN/SOYBEANS HISTOGRAMS 
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observations were compiled. Because the 8-bit conversion of the FSS data 
produces highly serrated histograms which disrupts visual inspection, the 
histograms presented in Figure 18 have been smoothed using a 3 count window 
moving average. 

The FSS histograms in Figure 18 and the NSOOl histograms in Figure 17 are 
quite similar in appearance. The mean and variance behavior of the corn and 
soybeans canopies display the same patterns in each data set. Note in 
particular, that the mid-infrared band distribution shows the same 
comparatively low within-crop variance when compared with the near-infrared 
band in the FSS data as was noted in the NSOOl data. 

Temporal Analysis of FSS Data 

The comparison of the FSS and NSOOl 8/30/79 data indicates that they 
provide common measures of the reflectance behavior of corn and soybeans 
canopies. The 8/30/79 data suggests that TM mid-IR observations may improve 
corn/soybeans discrimination on this data. However, analysis of current 
Landsat data shows that at this time of the growing season, after corn has 
tasselled, corn and soybeans are not difficult to separate using the visible 
and near-IR bands. 

FSS single scan observations for 6/11, 6/29, 7/16/79 were also processed 
to construct histograms of the corn and soybeans reflective patterns. Figures 
19-22 present the histograms of corn and soybeans FSS observations on these 
dates and the 8/30/79 observations. Table 9 provides the percent ground cover 
observed by the field enumerators on these dates. 

The histograms reveal a particularly interesting pattern of mid-infrared 
relectance for corn and soybeans when compared to near-infrared reflectance^ 

On 6/11/79 the FSS observes primarily bare soil and there is little difference 
in the comparative reflectance of corn and soybean fields in any of the bands. 



FIGURE 18 

8/30/79 FSS CORN/SOYBEANS HISTOGRAMS 
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TABLE 9 


PERCENT GROUND CO'VER FOR CORN 
AND SOYBEAN FIELDS IN 1979 


PERCENT 

COVER 

0 - 10 
11 - 20 
21 - 30 
31 - 40 
41 - 50 
51 - 60 
61 - 70 
71 - 80 
8,1 - 90 
91 - 100 


CORN 

SOY 

6/11 6/29 7/16 8/30 

6/11 6/29 7/16 8/30 

17 

30 

21 

6 

1 2 

7 

6 

10 1 

13 

2 3 1 

9 

1 3 

5 

1 6 


12 


7 

36 38 

3 33 
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CROP HISTOGRAMS FOR TM SIMULATED FROM FSS 
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CROP HISTOGRAMS FOR TM SIMULATED FROM FSS 
AgRISTARS ITS WEBSTER CO. IOWA 1979 CROPPING SEASON 

TM BAND 4 (0.76 -O.SOi&fn) 
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CROP HISTOGRAMS FOR TM SIMULATED FROM FSS 
AgRISTARS ITS WEBSTER CO. IOWA 1979 CROPPING SEASON 
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Conclusions 

The GlSS/Columbia Earth Resources group has been intensively involved, 
during 1981, with activities directed at simulating TM observations and 
evaluating the potential information content of these observations. The 
NASA/AgRISTARS field measurements program continues to provide extensive data 
to conduct this research. The GISS/Columbia staff have developed processing 
and analysis techniques to utilize this data base. The 1981 research program 
has successfully addressed several outstanding NASA research objectives 
concerned with preparation for the analysis of satellite acquired TM 
observations. The more extensive field observations collected at the 
AgRISTARS supersites during the 1980 and 1981 growing season present a major 
analysis challenge. The experience gained by the GISS/ Columbia staff in 
analysis if the 1979 AgRISTARS data during this year will be invaluable to 
exploitation of the more recent field measurements and in analysis of the 
soon-to-be acquired satellite TM observations. 

C. Analysis of Snow Cover and Wetlands from Landsat Data 

The GISS/Columbia group,, in cooperative research program with U.S. Army 
Corps of Engineers applications scientists at the Cold Regions Research 
Laboratory, Hanover, New Hampshire, have, during 1981, carried out analysis of 
selected Landsat scenes. This research objective is to determine the utility 
of Landsat numerical observations for mapping and analysis of snow cover and 
wetlands. Michael Miller, a graduate student of the Department of Geography, 
has served as senior analyst for these activities. Drs. Goward and Ungar of 
the GISS/Columbia group provide research direction in consultation with CRREL 
scientists, Carolyn Merry, Dr. Jerry Brown and Dr. Ike McKim. The following 
reports describe the research activities. 
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By 6/29/79, corn fields average 50% ground cover and soybeans 35% ground 
cover. In the mid-infrared corn fields are less reflective than the soybeans 
fields yyhereas in the near-IR the reverse is true. On 7/16/79 corn canopies 
are reported to have 100% ground cover whereas the soybeans fields average 65% 
but are highly variable. The mid-IR data show that corn is still less 
rel active than soybeans and better separated than on 6/29/79. However, in the 
near-IR the soybeans reflectance is highly variable and mixed with the corn 
reflectance. By 8/30/79 both crops are at 100% ground cover and soybeans are 
more reflective than corn in both TM bands. These histograms suggest that the 
TM mid-IR observations will significantly improve early season corn/soybeans 
discrimination when compared to currant Landsat observations. 

Discussion 

This investigation of the mid-IR reflectance behavior of corn and soybeans 
canopies suggests that TM observation in this wavelength interval will improve 
both early season and through-the-season discrimination of corn and soybeans. 
These conclusions are drawn, however, from one growing season in one location. 
Further analysis of observational data will be required to establish the 
general applicability of these findings. 

In depth evaluation of these results will require examination of the 
physical processes that lead to differential reflectance behavior of 
corn/soybeans canopies in the mid-IR. Interactions between leaf reflectance 
properties, canopy architecture and background soil need to be examined. 
Although further empirical studies may be able to confirm the behavior 
observed in this study, prediction of where and under what conditions this 
behavior will be observed will require physical understanding. 
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Progress Report: Remote Sensing of Snow 

Michael S. Miller July - December, 1980 

Introduction: Landsat im being investigated in conjunction with 

ground snow course measurements to determine the utility of the 
remotely sensed data for snow cover evaluation. Although 
satellite data proven useful for estimating areal extent 

of snow, applicability in forested hydrologic basins, and for other 
descriptors of snow cover (e.g. depth, water equivalent) are less 
obvious . 

Study Site: Allagash, St. John Rivers confluence area, Maine. Land 

Cover is a mix of deciduous, coniferous, and mixed forests, clearings, 
agriculture, and water surfaces. 

Achievements (July - Dec., 1980) 

1. Using a May 31, 19)Ef78 Landsat image, the 300X300 pixel study 

area (153 sq . miles) was classified using the GISS-MAPl classification 

routine. Results, in percent cover of the area^ are: 

Water 2 Mixed Forest 52. 

Deciduous Forest 16.3% Clearings . 5 % 

Coniferous 28.2% Agriculture .5% 

2 . Using these percents as weighting factors for ground snow depth 
measurements, average snow depths for the study region were determined 
for four days for which good quality (Landsat coverage is available; 


February 

11, 

1978 

29.01" 

January 

11, 

1979 

23.63" 

December 

24, 

1978 

20.94" 

January 

6, 

1978 

17.06" 


56 


original page 13 

OF POOR QUAUTY 


3. Landsat MSS digital data have been investigated to determine whether 
there is a crrelation between snow depth and the measured intensities 
in the four spectral bands. All sites were thus evaluated. Band 
ratios and summations were considered, as was the inclusion of a sun 
elevation factor to compensate for differences in incoming radiation 

on different days. In most cases (see Figure one a^d) little 
correlation between snow depth and landsat counts/energies was noted. 

4. Although radiances measured in reflected visible and near-infrared 
wavelengths within an individual pixel do not increase (or decrease) 
with increased snow depths beyond several inches, it is believed that 
the regional landscape does change. In the northeast forests, the 

U'' . 

vegetation canopy is d Mixtures of clearings and 
coniferous and deciduous forests at various stages of growth allow 
upper canopies and understories of differing heights and densities. 

As snow accumulates, different levels of vegetation would be expected 
to be covered with snow. Figure two a-d show band 7 grayscales for the 
300^300 pixel Allagash area for four days. As is seen, increasing 
snow depth results in increasing 'whiteness' of the area. 

5. Histograms of the energies have been investigated to quantify this 

change (figure 3 A no snow* scene in autumn shows a band 7 histogram 

that is unimodal and of low variance (3-a) .A?(s snow depth increases 
(3-b to 3-e) , there is a decreasing in peakedness and increasing of 
skewness toward the bright end of the distribution. This indicates 
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that as snow depth in a diverse land ccver increases, increasingly 
more pixels will become 'snow pixels', migrating in the regional 
histogram from the apparently snowless vegetation peak to the 
brighter (snow) intensities. 

Projected work (<Jan. ^ June 1981) 

Examination of gray scales and histograms for additional snow 
scenes of the ^»llagash/st . ^ohn area will continue. This is to 
dtermine the validity of the preliminary observations. 

In addition, gray scales will be examined in conjunction with 
the derived classification map to determine which vegetation types 
contribute to increased brightness as snow depth increases . 

similar examination will be initiated in the Danville, Vermont 
area in the Sleepers River Basin. Records of snow depth for a 
variety of land cover types are available for this intensively monitored 
basin in northeastern Vermont. 

Problems encountered to date: The major obstacle to this research into 

snow cover is the low number of available Landsat scenes for this region 
during snow seasons. This is due primarily to the satellite's cycle 
and to the frequency of cloud cover, particularly in the important 
snow melt period. 
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D. February 11, 1978 


Average snow depth: 29.01 


Figure two, continued. 


C. January 11, 1979 


Average enow depth: 23.63 
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Figure two. ^andsat HSS 

Band 1 , gray acalaa. 


A. January 6, 1976 

Average anow depth: 17.06" 


B. Decerr±:er 24, 1978 

Average snow depth: 20.94" 
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Figure Three. HISTOGRAMS: LANDSAT MSS BAND 7 
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ABSTRACT 


Each spring the SSARR (Streamflpw Synthesis and Reservoir Regulation) 
model is used iv> the Upper Saint John River Basin, Maine, to forecast runoff 
due to snowmelt and precipitation. One of the main Inputs to the model is 
the extent and distribution of the snowpack. Using point snow course measure- 
ments to determine the extent and distribution of the snowpack is a commonly 
recognized problem. A technique for using Landsat digital data to map the 
distribution of snow cover/land use categories has been evaluated for use in 
the SSARR model to improve the accuracy of spring runoff predictions. 

Landsat MSS digital data has been investigated to determine whether there 
is a correlation between snow depth and the measured intensities in the four 
spectral bands. Band ratios and summations of the four bands were considered, 
as was the inclusion of a uun elevation factor to account for differences in 
incoming radiation on different days, A direct correlation between snow depth 
and Landsat counts/energies does not occur. 

Although the Landsat radiances measured in the reflected visible and near 
infrared wavelengths for a single pixel do not increase (or decrease) uniforaly 
with increased snow depth beyond several inches, it is believed that Landsat 
radiance values on a regional scale do change. Histograms of the Landsat radiance 
values for a 300 x 300 pixel (153 sq. tni.) study area have been Investigated to 
quantify the change. A histogram for a no snow scene in the late autumn (27 October 
1978) exhibits a <;nlmodal distribution with low variance. As the snow depth 
Increases over the total landscape, there is a decrease in peakedness and a 
corresponding Increase of skewness toward the bright end of the Landsat distri- 
bution of radiance. This fact indicates that as snow depth in a landscape 
increases, additional pixels become 'snow pixels', migrating in the regional 
histogram from the apparently snowless vegetation peak to the brighter 'snow' 
intensities. 

Various indices of skewness and peakedness (kurtosls) have been examined to 
quantify the above observation. An index of potential use is defined asi: 


« 83 * pO-E- C tLo* ; . ? ! ) , *C0S(i) C 

^ pop (mode) 


64 


ORIGINAL PAM 18 
OF POOR QOAUTY 


in which !■ the index of ■nowneiSi *3 ii the aonent coefficnt of ekewneet, 
pop (mode) ie the pixel population of the histogram mode, pop (total) ia the 
total pixel population and COS(ft)eiB the cosine of the solar xeiii^ angle. 
Through linear regression techniques » It was determined for tiCe'sbo x 300 pixel 
Landsat area that: 


Snow Depth - 2.079*(lg) -1.111 

in which snow depth is measured In inches. The index of anowness, was 
found to predict mean snow depth with a standard error of estimate of 4 4.2 
inches at a 95% confidence level. For different mixtures of canopy cover* 
the slope and y-lntercept would vary. 

As a preliminary test of this index of snowness^two 

winter Land sat passes from previous years were examined for 

the Allagash - St. John region. Using the index and equation 

of regression^ mean snow depth for February 11* 1973 and 

April 19, 1974 were estimated at, respectively, 35.49*’ and 

o. 

15.69", - 4 . 44 ". During those years, ground snow course measurements 

were not available for the Allagash stations. Closest measurements 

were taken at Connors, New Brunswick, 12 miles east on the 

St. John River. Snow depths at Connors were: 

February 1, 1973 33.3" 

February 19, 1973 37.9" 

April 16, 1974 23.7" 

May 7, 1974 11.3" 

Xl. Bibliographic search for literature concerning the 
remote sensing of snow pack characteristics has continued 
during the report period. 


III. Visit was made on March 25 to the Sleeper's River 
Research Watershed. Snow courses, meteorological stations and 

streamflow measurement sites were examined in anticipation of 
studying the land cover - snow pack relationship for this basin 
in Vermont « Search of Lands? t scenes is underway. 


PROGRESS REPORT 


- Mlchasl 6. Mlllsr 

Project: Simeon, Nebraska - Wetlands Classification 
£or__Perio(3 ; September/October, 1980 

fibjegbiv^t To classify Landsat MSS digital data for June 12 , 1978 
(scena-i.d. 3099-16494) for Simeon (Nebr.) SW and SE Quadrangles. 
This area lies within the Valentine National Wildlife Refuge. 
Three cover types are classified: open water, marsh, and 
subirrigated meadow. Unclassified areas are 'classed* as 
uplands. This classification is performed for the Waterways 
Experiment Station, Army Corps of Engineers, 

Procedure: 1. The MAPI classification algorithm developed by 

Dr. S. Ungar of the NASA Goddard Institute for Space studies 
is used in this project. A description of this algorithm is 
provided in the appendix to this report. 

In MAPI, the digitised Landsat counts are converted to 
energy values (mW cm“^ sr“^) , For each classified cover type 
a signature, or group of signatures is required. The area 
within the Landsat multispectral space surrounding each signature 
and classified to its cover type is determined by user specified 
parameters: and WBRT (maximum angular delta, and brightness 

weighting) . 

2. Signatures were selected by the evaluation of test sites 


$ I 


within the larger study area. Six to eight test eitee were 
■elected from the quadrangles for each of the three cover types. 

This was accomplished by using the wetlands map overlay (scale 
1(24,000) provided by the Omaha District, Corps of Engineers. 

The overlay was superimposed onto a NSS band 7 grayscale printout. 
Test sites were delineated for each cover type using the ortho** 
photoquad for reference. The test sites ranged from 5 hectares 
(12 acres) to 40 hectares (100 acres) . To maximize 'purity' 
of cover type, the sites were selected from the interior of 
larger homogeneous areas. 

3. Unsupervised classifications were run to distinguish distinct 
classes within each test site.. Printouts of the Landsat (^rdf^rgies 
were used to derive signatures for each class. From the unsupervised 
classifications it was determined that two signatures were 
necessary to most accurately describe each of the cover types. 


Cover IVoe 


Sianature 

( in enerav 

; mW cm” 

2 ) 



Band 4 

1^-5 

ind__6^ _ 

. £and_7^ 

1. Water 

(a) 

.29 

.15 

.08 

.09 


(b) 

.42 

.25 

.20 

.21 

2 . Marsh 

(a) 

.41 

.26 

.37 

1.17 


(b) 

.47 

.32 

.50 

1.45 

3. Subirrigated 

(a) 

.47 

.32 

.58 

1.7i. 

Meadow 

(b) 

.46 

.31 

.61 

1.8^ 
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4. The signatures derived in the unsupervised runs were 


used in the supervised procedure for a larger area approximati.ng 
the Simeon SW Quadrangle. The supervised classifications were 
run with a variety of/Jj^^ / WBRT combinations, thus altering 
the set of pixels classed within each cover type. Considerable 
visual agreement between the classification and wetland map was 
achieved after fifteen runs. Refinements in agreement resulted 
after ten additional runs in which signatures were adjusted and 
rearrangements made in the order in which the signatures were 
assessed for each pixel (Note: order is of importance in MAPI, 
as the algorithm is deterministic rather than probablistic) . 

It was found possible to eliminate the need for one of the water 
signatures without a reduction in classification accuracy. 

5. In conversation (mid-September) with Horton Struve of the 
Waterways Experiment Station, it was decided to extend the 
classification to include the Simeon SE Quadrangle. Land cover 
within the two Quadrangles is very similar. Nevertheless, 
several additional runs were required to take into account some 
spectral variations across the region. 

6. Final signatures, in order of their evaluation within 
MAPI are: 
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Cover ‘Tvoe 

Siqnatur 

e 

» V 

WBRT 



^ 1 ^ 


- 2 ”” 


1. Marsh 

.47 

.32 

.50 

1.38 .041 

.1 

2. Sub. Meadow. 47 

.32 

.58 

1.71 .08 

.1 

3 . Water 

.42 

.25 

.20 

.21 .52 

.4 

4. Sub. Meadow. 40 

.30 

.58 

1.80 .12 

.1 

5 . Marsh 

.40 

.32 

.30 

.68 .12 

.1 


7. Percent cover was evaluated for each of the following 
classes: open water, wetlands (marsh + subirrigated meadow), 
and upJ, ».'nds for the entire Simeon SE, SW area. In the table 
below, these are compared to the corresponding figures for the 
Simeon SW area used at the Waterways Experiment station for a 
ground truth standard. Comparable figures for the entire area 
were not available. 

Cover Type Gr o und-Truth MSS Difference 


Water 

8.8% 

8.0% 

.8 

Wetland 

35 . 1 % 

32.6% 

3.1 

Marsh 


4.6% 


Subirrigated Meadow 


28.0% 


Upland 

55.5% 

59.4% ; 

3.9 


8. JA tape of the digitized final classification is being 
prepared for Horton Struve, Waterways Experiment station. 
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The following description of the GISS MAPI classification 
^ algorithm is taken from; 

Merry, C.J. et al., Preliminary analysis of water 
equivalent/snow characteristics using Landsat digital 
processing techniques ^ 1977, Cold Regions Research and 
Engineering Laboratory, Hanover, NH. 


The LAIIDSAT I*ISS ohservation nay be thought c-f as a point in a four- 
dimensional "color" space, where the values along each axis represent 
the radiant energy received by the satellite in one of the four bai.is 
(this is illustrated for three bands in Fig. 2a). Observations lying in 
a similar direction from the origin in this four-dimensional color space 
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are said ,to te similar in color regardless of their total radiant 
energy. The distance of an observation from the origin is a measvire of 
the total radiance associated with that data point. The algorithm is 
primarily designed to combine observations similar in color into the 
same classification category. There is provision for evaluating brightness 
differences and for weighting these differences in with the color dis- 
crimination when constructing the classification categories. 

Discrimination based solely on color is obtained when one examines 
the difference in direction between the color vectors (observations). 
If^the_angle between the observations is smaller than some user-defined 
Qvitevion the vectors are considered to be lying in the same direction 
and, therefore, the observations are placed in the same category. 

There are two modes in which this classification scheme may be 
employed, supervised and unsupervised. In the supervised mode the user 
specifies a signatm'e (the energy distribution in the four LAITOSAT bands). 
If an observation lies within an angle smaller than the user-defined 
criterion, it is said to belong to the category represented by the 

signature (Fig. 2b). Therefore, all vectors lying within a cone of 
angle, dnax* about the signature representing category X belong to 
category X. 

In the unsupervised mode the color vector corresponding to the first 
observation is compared with all subsequent observations. If color vector 
1 is similar in direction to color vector 2 (i.e, 66 £ observation 2 

is placed in the same category as the first observation (Fig. 2c). In a 
similar fashion observations subsequent to observation 2 are compared to 
the second observation and so on right up to the last observation. If in 
the process of constructing categories, a member is found which belongs 
to a previous category, the new category is chained to the original 
classification category forming one Joint category. In effect the un- 
supervised classification will form several categories based on a criterion 
specifying maximum color difference permissible bexiween members of the 
sane category. 

In addition tc discrimination based solely on color, the GISS al- 
goj-ithm provides the capability of weighting total radiance {alhedc) dif- 
ferences into the discriminant. The percent albedo difference between 
two observations is computed. This normalized albedo difference is then 
combined with the color difference angle (expressed in radians) by per- 
forming a weighted average in the Rj*lS (root mean square) sense. This 
albedo-weig ,ted quantity is now compared with the user-defined criterion. 

A relatlvei'^ small albedo weighting allows very large albedo differences 
to disqualify observations that are similar in color from membership in 
the same category, thereby adding a second level of discrimination. 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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Band 3 


ORIGINAL PAGE IS 

OF POOR QUALITY 



a. A color vector vhich is illustrated for three har.is. 


bnd 3 



h. Gutervised nicde . The user-defined criterion, defines category X about the 

specified sifrnature, §, Any color vector that lies within this cone belonrs to 
category X. 


■•nd 3 



'inourer vised tnode^ is similar in direction to (60 <_ and placed in category 

1. 0^4 is similar in direction to S 2 ^iid plared in category 2. However, "Si* is also 

similar in direction to B 3 (category 1). Therefore, category 1 is merged with category 2. 

Fig'ure 2. Concept of the four-dinensional ''color" space, represented in three bands. 
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Project; Burlington, Vermont Quadrangle: Wetland Mapping 

Objective: To map potential and definite wetlands; evaluate 

seasonal images. 

Procedure: Two Landsat images of the Burlington area were 

classified using the GISS MAPI classification package 
The dates of the two passes were 13 April 1977 (SCNID 
20812-14372) and 26 June 1975 (SCNID 5068-14433). 

The first stage in the classification of the two dates was 
the development of a standard routine for mapping potential 
wetlands. Potential wetlands are here defined as that set of 
p,i?iels in the scene which include all wetlands. That is, errors 
of omission are minimized. However, there may be large errors 
of comission. The second stage was to perform a second classifi- 
cation on these potential wetlands in which only definite areas 
of wetlands are identified; errors of comission are minimized. 

In this second stage, areas of open water, emergent vegetation, 
and forest/shrub-scrub wetlands were segregated. The two classi- 
fications used in conjunction may give a good indication on the 
locations of the wetlands. 

13 April 1977 

1) After examining numerous classification possibilities for 
potential wetlands, the approach finally selected was a one 
signature, two-band classification using MSS-4 and -5. In the 
unsupervj. jed mode of MAPI, pixels corresponding to open water 
were found to cluster around a signature of 0.44, 0.25, 0.00, 

0.00 when weighting of bands equalled 1.0, 1.0, 0.0, 0.0. This 
signature was then applied to the supervised mode with a range of 
increasing delmax values. The final namelist parameters selected 
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for this classification were SIG « 0.44| 0.225» 0.00# 0.00} 

DELHAX « 0.215} WBRT « 0.5f WBND « 1.0, 1.0, 0.0, 0.0. 

2) To classify definite wetlands three classes were use^l.! open 

water, emergent vegetation and forest/shrub-scrub wetlands. All 
areas identified as potential wetlands in step 1 were classified 
using a one-band (MSS-7) classification. The namelist paraimeters 
used were: a) water: SIG « 0.00, 0.00, 0.00, 0.17; DELMAX « 0.7; 

WBRT K 1.0; WBND « 0.0, 0.0, 0.0, 1.0; b) emergents: SIG * 

0.00, 0.00, 0.00, 0.53; DELMAX 0.25; WBRT « 1.0; WBND « 0.0, 0.0, 
0.0, 1.0; c) forest-shrub-scrub: SIG « 0.00, 0.00, 0.00, 0.55; 

DELMAX « 0.444. 

26 J;ine 1975 

3) Using the seune procedure as in (1), MSS-4, -5 signature of 
water (0.40, 0.15, 0.00, 0.00) was applied in the supervised mode 
with DELMAX = 0.3, WBRT -= 0.5, and WBND = 1.0, 1.0, 0.0, 0.0. 

4) The nature of the landscape in June made it possible to 

successfully classify only open water of the three wetlemd 
classes. In these cases a four-band classification most accurately 
identified water: SIG - 0.40, 0.15, 0.10, 0.11; DELMAX - 0.4; 

WBRT « 0.1; WBND * 1.0, 1.0, 1.0, 1.0. 

Comments 

1) Seasonal variations are highlighted here. Spring images are 
best for identifying wetlemds in the Northeast - water table 
heights are highest in this season. Additionally, deciduous trees 
have not leafed-out, allowing surface conditions to be viewed. 


ORIGINAL PAGE S3 
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During the summer, canopy closure is more con^plete and the soil 
is mostly dry. 

2) The class of forest/shrub-scrub wetlands in the classification 
includes deciduous species. Wet evergreen woodlands cannot be 
separated from dry evergreen woodlands in either season. 

3) The category of intermittent exposed streams as included in the 
ground truth map of wetlands cannot be identified with Landsat 
because of the features' small size. 

The Tape 

The tape is 9-track, 800 and contains 5 files 

separated by End-of-File marks. All records are 240 bytes in 
length. The first file contains a single "header" record in 
EBCDIC format which describes the subsequent files' characteristics. 
Files 2, 3, 4, and 5 contain the classification results stored as 
EBCDIC symbols. Each file contains 224 records. File 2 contains 
potential wetlands classification for April 13, 1977. File 3 
contains definite wetlands classification for April 13, 1977. 

File 4 contains potential wetlands classification for June 26, 

1975. File 5 contains definite wetlands classification on June 26, 
1975. 

For the potential wetlzmds classifications (files 2 and 4) 
the symbol '1' represents potential wetlands; '-' represents non- 
potential wetlands. 

For the classification of definite wetlands on April 13 (file 
3), '-' represents non-wetland, '!' represents open water, '2' 
represents emergents, and '3' represents fore st/shrub- scrub wetlands. 

For definite wetlands on June 26 (file 5) , only open water 
is identified as *1'. Non-wetlands are represented by '-'. 
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Conclusions 

These investigations have contributed to the development of U.S. Army 
Corps of Engineers applications of current Landsat data. Highlights of the studies 
include the GISS algorithm performing the best is a comparative study of 
classification algorithms on the Nebraska Wetlands study and the innovation of 
a new analysis concept for snow cover applications. This cooperative effort 
has proven highly productive and will continue during 1982 with emphasis being 
placed on refinement of the snow cover analysis technqiues. 

D. Agro-Environmental Application of Satellite Remote Sensor Data 

During 1981 GlSS/Colunibia staff members continued pursuit of additional 
applications of satellite remote sensor data to selected projects of personal 
interest and that support their educational and professional objectives. In 
general, these studies address specific questions of interest to geographers 
which will also contribute to increased knowledge of numerical remote sensor 
data. Presentation of this research has occurred at regional and national 
meetings of the Association of American Geographers as well as at the 
NASA/NOAA cosponsered Remote Sensing Educator’s Conference (CORSE-81) held at 
Purdue University (Appendix A). 

Landscape "Greenness" 

Helene Wilson has proposed that Landsat observations provide a unique and 
universal measure of landscape vegetated condition, called landscape 
"greenness". Considerable research has been conducted in agricultural 
applications of Landsat data that shows a strong relation between various 
attributes of crop canopies (e.g. , percent ground cover, leaf area index, 
biomass) and various transforms of the Landsat visible and near-infrared 
observations. Although it has not been proven, this approach to Landsat data 

would appear to be generally applicable to all landscapes. That is, the 
degree to which a landscape is vegetated should be related to its "greenness" 
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as observed by one or more vegetation Indices derived from Undsat 
observations. 

Helene has selected a study site In the Hartford, Conn, region, where a 
complex variety of land cover conditions Including forest, agriculture, and 
urban are present. Her hypothesis Is that regardless of the cover type, one 
or more of the Landsat-derlved vegetation Indices will Indicate the degree to 
which the landscape Is occupied by photosynthesizing vegetation. If this Is 
true then Landsat data can be used as a universal measure of landscape 
conditions where the vegetation "greenness" may serve as an Indicator of 
current environmental conditions. 

During 1981 Helene acquired low altitude color Infrared stereo photography 
and conducted a field measurements program to evaluate the relation between 
ground conditions and the air photo observations. The aerial photography will 
provide Landsat subpixel measures of percent area covered with living 
vegetation and vegetation height. Through intensive analysis of the 
photography she will establish the vegetated conditions within Landsat 
observations of various land cover conditions. Through comparison of the 
photo-derived observations and the Landsat Data, she will test the research 
hypothesis. 

Albedo Studies 

Chris Brest, during 1979 and 1980, collected field measurements of the 
reflectance of building roof tops and parking lots in the Hartford, Conn, 
region. These observations are being used to calibrate a set of 28 Landsat 
observations of the Hartford region acquired between 1972 and 1978 and 
covering eleven months of the year. Calibration of these data to percent 
reflectance and the combination of Landsat visible and near-infrared 
measurements provides a good estimate of surface albedo. This further permits 
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analysis of urban-rural albedo differences as a function of season and cover 
type. Albedo has been suggested as a significant factor In causing 
urban/rural clltnate differences. Chris's analysis will evaluate this 
hypothesis. 

A preliminary analysis of the data (Appendix A) has confirmed that the 
calibration procedure is effective (error rate of 2 ~ 3 %) and that urban-rural 
albedo differences may be smaller (approximately 5 %) than Is frequently 
reported. This appears to be related to the Inverse variation of 
v1 si bl e/near-infrared reflectance through the season; a phenomenon generally 
not Investigated In albedo studies. During 1981 Chris completed analysis of 
the entire 23 Landsat observations and is currently drafting his research 
results for submission to the Department of Geography faculty as his doctoral 
dissertation. 

Michael Miller has approached derivation of surface albedos from Landsat 
data through numerical modelling of atmospheric radiative transfer. By 
calculating the effects of scattering, as described in the U.S. Air Force 
Cambridge Laboratories research, and absorption, as described In the 
Smithsonian Tables, he evaluates the relation between ground reflectance and 
observed Landsat radiance. Seasonal variations In atmospheric conditions are 
considered but not short term meteorological conditions. 

Michael has applied the model to Landsat observations of eastern Long 
Island, N.Y. where diverse land cover conditions present a variety of surface 
albedos. The results, although promising, suggested that his method for 
handling atmospheric absorption of radiation is less realistic than desirable. 
The albedos calculated are the right order of magnitude but appear to be 
overestimated. Using only the scattering model he derives more reasonable 
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figures. He is currently drafting his research results for submission, as h1s 
master's thesis, to the Oeparttnent of Geography. 

Kenyan Agricultural Systems 

Tina Cary is conducting a study of agricultural systems in western Kenya 
(east Africa). Her hypothesis is the Landsat observations record regional 
variations in the agricultural systems present. In particular, she proposes 
that, because of striking regional variations in the intensity of agricultural 
development and regional variations in and interaction between cultural and 
physical factors that effect local agricultural practices, Landsat 
observations may be used as an indicator of the agricultural system present in 
a location. 

To identify the relation between regional variations of agricultural 
systems and Landsat observations she will examine the numerical properties of 
spatially aggregated Landsat data. Through analysis of the data at various 
levels of aggregation she will seek to isolate those scales (or resolutions) 
at which various cultural and/or physical factors are most highly related to 
landscape appearance. Tina will use this information to develop a systematic 
approach to landscape regionalization which will permit analysis of the 
regional patterns of agricultural systeins. 

During 1982 Tina will travel to Kenya, under a Fullbright scholarship, to 
conduct field observations to support her analysis. She will be in the field 
for one year traversing the 10,000 square mile area that represents her study 
site. The study site consists of the nominal ground coverage of one Landsat 
scene and thus represents an ambitious field project. Because of her previous 
field experience in Kenya with Dr. Frank Conant of Hunter College, Tina 
should be highly successful in accomplishing her objectives. 
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Conclusions 

The diverse research projects undertaken by Columbia students in 
agro-environmental research provide major contributions to the GISS/Columbia 
research program with a minimum cost to NASAt The activities in general* do 
not draw heavily of the group resources and often, as in Tina's case, are 
significantly supported by other funding sources. However, the outcome of the 
research significantly contributes to advancing understanding of satellite 
earth observations. The success of these projects are dependent on the 
working relation between GISS and Columbia scientists and demonstrate the 
productive nature of the cooperative agreement. 
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DyNAMICS OF SURFACE ALBEDO by Christopher L. Brest 
’his study assesses the role of aTbedo as a contributing factor 
5 the formation of urban heat Islands. Field spectroradlometric 
measurements are used tc calibrate tMenty>sevcn Landsat observa- 
tions of the Hartford, Connecticut region. Calibration of these 
remotely sensed data permits analysis of the spatial and temporal 
variations of surface albedo for the entire region. The study In- 
vestigates albedo differences between man-made urban surfaces and 
natural vegetative surfaces. Additionally, spectral reflectance 
In the visible and near-infrared are being studied to determine 
their contribution to net shortwave albedo. 

SEASONAL TRENDS IN LAND COVER ALBEDO ON EASTERN LONG ISLAND by 
Michael S. Miller 

Land cover albedo Is evaluated using multitemporal Landsat MSS 
data. This Is accomplished by applying the Elterman atmospheric 
attenuation model, Temporal variations in Incoming radiation at 
the surface are first determined with the model and compared with 
ground pyrenometric measurements. Assuming 100 percent reflec- 
tion, the model is used to estimate radiation reaching the satel- 
lite sensor. The ratios of Landsat energy measurements to these 
estimates are used to derive albedo values. The Riverhead, 

Long Island study site Includes agriculture, oak forest, pine 
barrens, water, and urbanized lands. These lend cover types are 
examined for spatial and seasonal variations In surface albedo. 


PATTERNS OF AGRICULTURE: A FUNCTION OF SCALE by T1na K. Cary 
Studies in agricultural geography have general ly' occupi ed two ex- 
tremes on a scale continuum: detailed field mapping, and large-area 
regionalization. The opportunity now exists to Investigate spatial 
patterns of agriculture at meso-scales as well, using the Landsat 
data base. Tn this study, patterns of agriculture In the tropical 
highlands of western Kenya are being investigated for a range of 
scales. Landsat digital data are spatially aggregated to generate 
the range of scales, Using this approach, the gap in the scale con 
tinuum can be bridged, permitting explicit consideration of how pat 
terns of agricultural land use appear to vary with scale. 


LANDSAT MEASUREMENT OF LANDSCAPE GREENNESS by Helene Wilson and 
Samuel N. Goward 

The principal Tnformatlon contained in Landsat mul t1 spectral data 
is the degree to which any landscape Is occupitd by photosynthe- 
slzlng surfaces and Is expressed In the relative differential be- 
tween visible and near-infrared radiances. Numerous studies have 
demonstrated the utility of vegetaiton indices based on this dif- 
ferential for assessing a variety of land cover conditions. This 
investigation extends previous work by explicitly relating the 
amount of actively producing vegetation cover within a pixel to r.? 
coroed radiances for a heterogeneous lands'tape. The "greenness" 
approach offers potential as a nomothetic solution to the probleift 
of exploiting the uniquely global features of Landsat observations 
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Paper presented at 1980 Annual Meeting, 

Association ot American Geographers Middle States Division, 
University of Delaware, Newark, Delawarf^, October, 1980 


OBSERVATION OF URBAN/RURAL ALBEDO CONTRASTS 
Christopher L. Brest and Samuel N. coward* 

Columbia University 

^Columbia University and NASA Goddard Institute for Space Studies 

INTRODUCTION 

Since albedo is a significtnt determinamt of net radiation it has 
been suggested that the modified urban climate may be due to the differences 
in albedo of urban and rural surfaces. This research is an investigation 
of variations in surface albedo in the Hartford, Connecticut region. The 
intent of the study is to assess urban/ruTal surface albedo differences 
and to ascertain the seasonal variability of these differences. 

It is hypothesized that the albedo of a mid-latitude city in a humid 
climate is significantly lower than that of its surrounding rural environs. 
The higher albedo of the rural area is attributed to vegetation in the 
summer and snow cover in the winter. It is further hypothesized tliat the 
urban/rural differences vary seasonally due to changes in rural albedo. 

The methodology employed in this study allows measurement of the 
albedo of the entire study site, not just 'representative' seunples. Urban 
areas are mosaics of various man-made and natural surfaces amd only by a 
study of the entire region can an understanding of the spatial amd tempoi'a! 
characteristics of albedo be acquired. Until recently equipment and tech- 
niques have not been available to undertake the type of study proposed here. 
Use of remote sensing techniques avoids the .problem of spot measurements 
or representative samples. The utility of remotely sensed, image formatted, 
data to measure solar and terrestrial radiation values has been demon- 
strated 
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AVAILABLE DATA 


The primary data for this study consists of 25 observations of the 
Hartford, Conn, area acquired between 1972 and 1976 by the Landsat multi- 
spectral scanner. The data are available in digital form and are being 
analyzed at NASA's Goddard Institute for Space Studies. The 25 observations 
give excellent seasonal coverage. The number of observations per month 
axe shown below: 

JAN. FEB. MAR. APR. MAY JUN. JUL. AUG. SEP. OCT. NOV. DEC. 

2 32 11324241- 

Additional data includes the radiometric measurements of selected 
urban targets collected by the authors which are to be used in the cali- 
bration procedure. 


RESEARCH DESIGN 

To successfully determine surface albedo from spectrally selective, 
remotely sensed satellite observations tliree questions must be addressed: 

1) How do you derive albedo from spectral measurements? 

2) How do you convert radiance to percent reflectance? 

3) How do you account for atmospheric effects? 

Albedo is defined as the ratio of reflected to total radiation in the 
short wave portion of the spectrum (<4.0 ym)®. It is therefore necessary 
to consider the reflectance of the surface in this entire range. Spectral 
intervals must be chosen which are representative of reflectance in this 
range. This necessitates a knowledge of the reflectance characteristics 
of the surfaces involved. Man-made surfaces such as concrete and asphalt 
display a fairly constant reflectance with changing wavelength. Vegetation 
has a very distinct reflectance curve, with a significant rise centered 



at about 0.725 vini separating the response of vegetation into two distinct 
segments: low reflectance in the visible and high reflectance in the near- 

infrared. This Variation in reflectance of a vegetative surface is accounted 
for in this study by Constructing a weighted average of the reflectance in 
tlie visible and near-infrared. The weighting factors, 60% for the visible 
and 40% for the near -infrared, are derived from the proportion of incident 
Bolaur radiation which lies below and above the vegetation rise at 0.725 ym®. 

The Landsat multispectal scanner operates in four spectral intervals; 
0,5-0. 6, 0,6-0. 7, 0.7-0. 8, 0. 8-1.1 ym. The use of two Landsat bejids which 
straddle the rise in vegetation reflectance can be employed to construct 
an accurate albedo measurement from Landsat spectral observations. 

The issues raised by the last two questions stated edjove are accounted 
for by the calibration procedure used in this study. The methodology employed 
utilizes a radiometer which provides multispectral reflectance readings 
of selected sites in the form of percent reflectance. The radiometer was 
designed specifically to provide these readings in support of remote sensing 
missions and is described elsewhere^. Once computed the percent reflectance 
of the targets is independent of atmospheric conditions. This is an impor- 
tant feature because once a target has been measured, it can be utilized 
in the calibration of saitellite observations from different passes. It 
eliminates the need to have ground based measurements taken coincidently 
with each satellite pass. A necessary condition for the success of this 
procedure is that the targets selected must be temporally stable. The 
targets chosen are large man-made surfaces: builclng roofs and parking lots. 

This will test a new type of target not previously used in calibrating 


satellite data. 



For this preliminary paper 14 targets whose reflectances were mea-^ 

•ured in the summer and fall of 1979 are being used. Field %a3rk is 
continuing in the summer of 1960 to increas^e the number of calibration 
targets. Th<i! 14 targets range from dark to bright; percent reflectances 
from approximately 4 to 50%. 

The radiometrically observed spectral reflectance meas^irements are 
used to calibrate the satellite data. Regression analysis is used to 
determine the relationship between ground and satellite observations of 
the targets sites (Fig, 1), Fach band is calibrated separately to transform 
the satellite observed radiances to reflectance and to remove atmospheric 
effects. The two bands are combined according to the weighted average 
previously discussed. Calibrated albedo maps can now be produced either 
in the form of computer printouts or video display. 

PRELIMINARY RESULTS 

Using the data from the 14 targets ■ r> began to evaluate our research 
design using 4 of the 25 satellite passes as representative of the seasons. 
The statistical relationships are shown below; 

TARGET LANDSAT STATISTICAL RELATIONSHIPS 
(Measured percent reflectance vs. Landsat radiance) 

MAY AUG NOV FEB (No snov'! cover) 


EXPLAINED 

VARIANCE 

BAND 4 
(Visible) 

0.90 

0.89 

0.90 

0.90 

(K^) 

BAND -6 
(Near IR) 

0.83 

0.86 

0.88 

0.86 

STANDARD 
ERROR OF 

BAND 4 

2.6 

2.8 

2.6 

2.6 

ESTIMATE band 6 

(% REFLECTANCE) 

3.9 

3.3 

3.2 

3.5 
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The values (0.B3-0.9) are highly encouraging and Indicate the 
potential of our procedure for calibrating satellite data. The standard 
error of estimate is higher than we had hoped but we arc confident that 
this can be lowered to approximately 1% through the addition of more 
calibration targets. 

Determination of the albedo of natural versus man-made surfaces was 
conducted using the USGS landuse-landcover classification 1:250,000 scale 
map of Hartford. For purposes of this paper we have selected six land- 
cover types to demonstrate our preliminary conclusions. They are: 

DR = Dense Residential A « Agriculture 

LR = Low Density Residential F »= Forest 

CD ® Commercial, Industrial, Transportation WN ■= Wetland, Nonforest 
(Downtown) 

In general, albedos of man-made surfaces are lower than those of 
vegetative surfaces. The seasonal variation of albedo for the six land- 
cover types is shown in Fig, 2. The variation is the greatest for the 
landcover types which have a high proportion of vegetation, i.e., A,F,WN, 
LR. The lemdcover types which are comprised mainly of man-made surfaces 
are lower and have much less seasonal variation. The driving force behind 
these albedo differences is the high infrared reflectamce from vegetation 
foliage. Reflectance in the near-infrared also shows seasonal variation. 
It is this high reflectance in the infraured which produces the summertime 
pea)? in the seasonal albedo shown in Fig. 2. This point is further rein- 
forced by looking at Fig. 3 which shows the ratio of BAND 6/BAND 4. There 
exists a high ratio in the spring and summer. The LR landcover type falls 
into this grouping because of the large proportion of vegetation present 


between the man-made structures. 

The three landcover types wMch might be considered 'urban' are 
shown in Fig. 4. This figure shows the individual band reflectances for 
the three types. Again the type with the highest proportion of vegetation « 
LR> has a much higher reflectance in the infrared during spring and suniner 
and is more variable. Notice the relative uniformity of the two landcover 
types comprised mainly of man-made surfaces, CD and DR, during the course 
of the year, 

Color video displays of the computer generated albedo maps clearly 
demonstrate the considerable seasonal dynamics of landcover albedos. 

Rural, vegetative, surfaces have a high albedo in spring and summer. 

Urban areas, or those comprised mainly of mem-made surfaces, have a lower 
albedo and remain more consistent throughout the four passes. In the fall 
and winter rural surfaces return to a lower albedo, similar in magnitude 
to that of urban areas. 

At the present time the albedo values derived are low in comparison 
with those stated in the literature. The reason for this may be the 
calibration panel which is used to derive the percent reflectance values 
for the targets. We believe the furnished percent reflectar»ce value for 
the panel is too low. At the conclusion of the s\immer 1980 field measure- 
ments program the panel will be sent out for laboratory testing to confirm 
this suspici,::n. 


CONCLUSIONS 

Results generated demonstrate that seasonal variation in rural albedo 


is significamt. To understemd urban/rural albedo contrasts seasonal 


variations must be studied. Certain technical limitations In the cali- 
bration procedure need to be resolved before the absolute stagnitude of 
urban/rural albedo contrasts can be assessed. The metliods applied in this 
research show great promise in providing appropriate data for analysis 
of both temporal and spatial variation of surface albedo. Final results 
of this research should provide a valuable contribution to studies of 
variations in interface energy budgets that result from variable land- 
cover conditions. 
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ABSTRACT 


In this study an exanination is made of spatial and temporal 
variations of urban/rural albedo in the Hartford, Conn« region. Such 
contrasts have been proposed as one cause of observed differences in 
climatic phenomena. Due to the complexity of ground cover, the deter- 
mination of spatial variability of surface albedo is difficult with 
with traditional ground based observations. In this study calibrated 
Landsat acquired remote sensor measurements are used for albedo assess- 
ment. Results indicate that urban/rural albedo differences vary 
temporally. This results primarily from seasonal changes in the 
reflectcince properties of vegetation canopies. Greatest albedo 
differences between urban and rural areas occur during the summer 
due to the high near-infrared reflectance by vegetation. The smallest 
differences occur during the winter, with the absence of snow. 
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THE problem 

Ceogrophy currently ploys o major role in remote sensing education 
in this country (Estes, et ol,, 1977), with the emphosls todoy on 
trolntng students in the application of remotely sensed obserwotionB to 
geographic problems, This emphosis reflects o prevailing view of 
remotely sensed doto i n te rp re to t i on os o well-established body of 
techniques to be opplled In solving o brood spectrum of geographic 
problems. This view Is also opporent in what we see os the current 
emphasis in geographic remote sensing research on applications studies. 

The underlying ossumptlon of this applications orlentotlon is that the 
informotion content v f remotely sensed data is known ond coincides with 
the types of informotion required for porticulor purposes. 

There is some question as to the degree to which onolysls of 
contemporary remotely sensed data for applications purposes is 
so t i B f oc to ry , The categories of Informotion typicotly of Interest in 
opplicotions work (e.g., Anderson, et ol., 1976) have not proved to be 
i n te rp re tob I e from current remotely sensed doto with conilstency ond 
confidence. Thie sltuotlon roises doubts pi to the degree to which the 
information contolned In the current gensrotion of sensor obserwotlons 
of the eorth's surface is understood. 

it 

SOURCE OF THE PROBLEM ^ k 

Modern remote senelng eystema hove been criticized Uy some y 

Investigators os not being effective In providing the types of I 

informotion for which ths use of troditionol cbservction eyitems (e.g., | 

oeriol photography) hos been proven successful. While this criticism I 

moy be viewed os merely o coneervotTve reoction to the oppeoronce of new I 

technology, it is o vo I t d criticism. The odvent of remote sensing from I 

spocs occurred during o period in which American geogrophicol | 

pho to t n terp re to ti on work emphoetzed opplicotions (Stone, 1 974). | 

Investigators oppeor to hove turned to the new doto oe another means of | 

solving the tome types of geographical problems to which they had been | 
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•pplying ovrfol photogrophy. Th«r« wok no phoo# of thworottcol 
Invootl flotl on of ihoto nww doto cowporobl# to thot which hod procodwd 
thw oppllcotlono phooo In oorlol photogrophlc tnlorprotatton . 

Not kurprlilngly, tnvootigotoro hovw oncounttrod d.fffcultUo fn 
•ddrwiotng o wot of probltnt thot hovt romotnod tho oomo with doto thot 
hovo ehongod liorkodty. Th# opotrol, toinporol, ond opoc tro I /rod t omo tr J c 
poremotoro of now oboorvotlon oyotami oro o 1 gn 1 f I con t ly dlfforont from 
theoo of whot moy new bo rogordod oo ” trod t t1 ono I " oyotoma. Tho mojor 
innovotiono ooooetotod with romoto oonotng from opoeo oro oo rpllowo; 

1) opotlol; uniform oboorvotlono oro oequlrod on o global ocolo. 
ond oro mort goneroilzod; 

25 tomporolJ covorogt to of o I gn I f I c on t I y higher frequency; 

3) opoc tro I /rod t omo tr i c : oboorvotlono oro numorlcol mooouromonto , 
In mony wovolongth bondo. 

Thtoo footuroB hovo o number of Impllcotlono for tho oKtroctlon of 
googrophtc Information from thooo now doto. For OKompIo, tho opotlol 
roioluttono typlcolly oncountorod In thooo now doto ofo cooroor then 
thooo choroctorlottc of moot oorlo! photogrophy. Tho rolotlonohfp 
botwoon ImogoB In photogrophy ond footuroo on tho ground io foirly well 
understood oo o rooult of occumulotod oxporlonco In OKtondlng vlouol 
perception of tho ourfoco to oorlol per opoc tl voo . However, there lo no 
rooion to oooume thot the ooooclotlono which con bo mode ot ono ocole 
ore volld ot onothor. Furthermore, oonoor oyotomo do not nocoooortly 
gonorolito londocopoo In woyo c onp lotont with our conceptual or 
eortogrophlc go no r o i i xa 1 1 on o . Slmllorly, the Incroooed tomporol 
rvoolutlon pooslfelo with ontolllto oboorvlng oyotomo fo quite dlfroront 
from thot of ” trod I tl ono I " oyotomo. Aerial photographic mloolono oro 
uoually plonnod for o time of yoor conoldorod t;> bo optlmol for o 
portlculor oppllcollon (o , g . , tree I do n 1 1 f t eo 1 1 on) , ond ropoot eovorogo 
lo rorely obtained more thon once every five ycoro. Ao o rooult, tho 
oolocilon of on opproprlote oboorvotlon dot® from omong tho many 
ovollobto for © given googrophlc locotlon becomoo prob S omo tlco I with 
thooe now data. 

AnolyolB of remotely oensed doto In numorlcol formot hoa been put 
forword os O^poaslblo oolutlon to oome of the difficulties presented by 
tho chorocteriotics of these new oboorvotlono. In portlculor, onolysls 
of the numorlcol doto ot the resolution limit of the sensor and 
utlltxotlon of Oil Wovolongth bonds slmultonoously has been expected to 
provide Information dbout the surfoce In greater detail thon thot which 
hum on beings ore oblo to oxtrqct from visual ossossmont of imagery. 
Automoted tnterpretotlon hos otso been seen as o moons of efficiently 
uitlixlng 0 doto bose thot hos been occumuloting with unprecedented 
rop I d I ty . 

Tho spectrol stgnoture porodlgm, widely employed In c omp u to r^boso d 
Interpretotlon of the multlvorlote, numorlcol doto, hos not mot with 
resounding^ success. For example. o survey (Jayroe, 1978) of 
( n v e a 1 1 g 0 1 1 0 n B using dlgltol Londoot doto to classify lend cover 
Indlcotes thot closslflcotlon occurocy figures ore on the overoge 
slgniflcontiy lower thon the minimum criterion of 86 percent often used 
for vlsuoi tmoge I n te r p r e to t I o n (Anderson, et ol,, 1976), Furthermore, 
spectrol signoture methods give occeptable results only over very 
limited oreoB ond only ot certoln times. In oddltlon, these methods 
noorly olwoys require slgnlflcont omounts of onclllory ground 
Informotlon, As o result, automated 1 n to r p r e to 1 1 on has yet to be proven 
os a meons of exploiting the globol ond temporol coverage offorded by 
sotelllte doto. Apporently, thw numerlcol onolysls techniques In wide 
use ore unobte to oddress effectively the probleme rolsed by the 
feoturee of contemporory remotely seneed eorth ob se rvo 1 1 one , 

We believe thot others, like ourselves, hove found that the results 
of onolysls of these new doto hove not conformed to Inltlol 
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•Kpvetotlonii How«v«r, «• do not b«lf«v« thoi th«i« now obitrvoifon 
•yotomi ohoutd bw rojociod, if fuggooi thoi' tho problom IlfO In our 
llwiifd undfratond I ng of tha netura af tha doto naw baing aequirad. A 
aetnprahanal va undaratohdlng af tha apactral, apetloj^ or tamporal 
ftroparilaa of londaeopaa, aa abaarvad with eentanporory aorth 
aboarvotlon ayatama, hoa not yat baan davalopad. 

iHAT NEEDS TO BE DONE? 

• Wa aaa o naad for boale raaaoreh to datarmtna tha Inharant 
Informotlon contant of tha doto. An opprooch that aa propoaa la to 
conaldar ranotaty aanaad doto oa o maoaura of ena or mora unknown or 
ioor ty-undaratood londacopa ettrtbutaa. Tha hypothaaai for thia 
raaaoreh fhoutd orlginota In OKomlnotion of tha romotaly nanaad doto 
rothar thon from aKtarnol objactfvaa, Tha gool of thio rtaaorch la to 
aaiobilah modala of londacopaa that may ba uofd to axplotn tha 
Informotlon contant of ramotaly aanaad data. 

Ona Hypothaala thot ahould be Invaatlgotad la whathar ramotaly 
aanaad obaervotlona provide maoauraa of contlnuoualy dIatrTbutad lond* 
ocopa phyalcol ottrlbutaa -- for axompla* molatura content, thermal 
Inertio, ond chioroynyll contant rother thon atmply Indicoilng the 
preaenca end eKtent of dlacrate cotegorlea or objecta. To teat th t a 
hypotheala requirea o thorough knowledge of the Interoctlona between 
a I ec tromogne 1 1 c energy ond londacopaa, oe eKpreaaad In the 
apec tr 0 I /rod I ome tr I c component of remotely aanaad data. Study of theaa 
Interoctlona la prevalent In remote aenalng reaeorch conducted In the 
npn-’photogrQphle regions of tha electromagnetic apectrum. For akomple, 
microwove remote aenalng reaeareh la concerned with the Influence of 
molaiure on obaervotlona: thermal propertiea ore conaldered to be key 
londacape ottrlbutaa In awploTnlng ihfrmal Infrared obit rvq t I ona . ThIa 
opprpoch la not oftan tokan In tha pack energy ronge of the aolor 
apectrum (vlalble ond neor> Infrared), becouae doto auch oa ihoae 
acquired by the Londaot ayatem iand to be viewed oa meraty on extonaion 
of whot la Already known from experience with photogrophlc obaervotlona. 
However, we believe thot thia preaumpilon la not correct. New remotely 
acnaed doto connot be viewed aimply oa o direct eutenalon of our vlauol 
eopobt l ltteo. 

Addltlonol atudy la needed to tnveatlgote the relot ‘onahlpe between 
!<*«■ rod'ometrlc measurement of londacope ottrtbutea ond spot I cl and 
t<iimporol foctora. For axompla, ona hypotheala might be thot the obtltty 
to meosure o porttculor londacope ottrlbute is Invoriont os o function 
of the epotlol resolution of the observing ayatem, Although one might 
•Kpect different ottrlbutes to be meoauroble ot different acoles of 
obaervotlOn, relotivciy little la known obout the woys in which aenaora 
generollze the aurfoce os o function of obaervoilon cell size. If the 
ability to meosure a portlculor attribute Is found to vary with aeolc, 
then Q new focua for I nvea 1 1 go 1 1 on becomes the de term I no 1 1 on of the 
aeoic ronge over which the ottrlbute ia meoauroble. 

With respect to temporol foctora, on Inltiol hypothesis that ahould 
be eKomlned la thot o given londacope ottrlbute is meoauroble ot oil 
times. The conalatency ond precialon of thot meosure may In foct vory 
with time -- for eKompIt, becouse of voriutlons in foe tors eKternol to 
the londacope, such oa sco^.4onol d'fffcrsnces In Intensity of aolor 
rodiotion. 

WHO SHOULD DO THIS RESEARCH? 

Speciollats In many disciplines hove contributed to odvonees In 
remote scnalng of the earth's aurfoce, Including engineers, 
mo theme 1 1 c T ena , phyalciats, ogronomiats, geologiata, ond othera. 
However^ i n vc a t i go to ra in eoch of ih-aie disciplinea undaratond ond ore 
eon earned with only o port of the observed londacope end/or of the doto. 
The deto, in controat, ore hollattc. In thot they represent Tntegroted 
obaervotlona of londacopaa. Geography la the dlaclpllne that clolma on 
Integroted opprooch to londacopea (Fennemon, 19)9), Therefore, 
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Qflog<-Qphcr6 beor o major ro »p on» I b II I ty ' for boBic r«ooorch on the 
romotoly sonsed obiervotion of londscapoB. 

WhMo .jr contontTon obogt the untquo opprooch of googrophorn would 
meet ttle d I aog ro«mo n t , the conoluoton we drow with reapect to 
the r • geogrophora anould be toking In boalc remote aonalng roaeorch 
la parlitis'l. leaa widely occopted. Ceogrophy oa o diaclpline oppeora to 
view It««lf largely oa o paaaive uaer of remotely aenaed doto, 
prlncipolly Intereated In ualng the doto oa on odjunct to current 
roBoarqh objectivea, The development of onolyaie techniquea, theorlea, 
ond rfpec I f too t i one for future ayatom deeigna fa left to other 
diaciplines. Mony geogrophere view boaic remote aenaing reaeorch oa 
"non-geogroph i c „ “ aince It to aeen oa being concerned with "technique" 
rother thon with whot la regarded oa aubatontivo geogrophic inquiry. 

This ottttudo hoB hindered geographers In moking remote aonalng 
reaeorch contributions commenaurote with the breodth of geogrophera’ 
perspectives on londecopaa. Wo cannot depend on othora to develop 
theories ond porodigma for ua. The concerna of the ayatemotlc diaci- 
plines ore not nocossarlly coincident with our own, Ceogrophy la the 
only discipline which takes on Tntegroted approach to the explonotion of 
oreol d i f f e re n t i 0 t i on on the eorth'a aurfoee. We routinely apply this 
opprooch in ottempting to exploln field obaervotlons ond ground 
measurements. Remote aenaing ayatems provide new observations ond 
mcosurements of oren' d ! f f o re n 1 1 o t i o n . We must assume our 
responsibility to seek .lew theorlea to explain the doto ond to provide 
others with the benefits of our insights. 


SIGNIFICANCE OF THIS APPROACH^ 

Remotely aenaed ofaaervations of the eorth'a aurfoee rolae quest’oris 
thot ore porticulorly geogrophic in their form ond scope. Geographers 
hove contributed a i g n i f I c o n t I y to development of ooriol 
pho to i n te r p re to t i on techniques by detoiled exominotion of the 
information contoined in the photogrophs. The new generation of 
remotely sensed observations hos yet to be exomined with such rigor. 
This foilure to conduct the needed basic reaeorch hoa constrained our 
ob’lity to extract and opply the geographic informotion rontoined in 
these new doto. We contend that the reaeorch npprooch we ore proposing 
will not only improve opplicotions of these new observotions to 
geogrophic problems of current interest but will o I so provide new woys 
of exomining londscopea that moy directly contribute to odvoncement of 
geographic methodology and theory. 


The d e m 0 n a t r 0 1 1 on thot remotely sensed observotions of the eorth'a 
aurfoee consistently measure selected attributes of (ondscopes should 
enobie the construction of models thot describe the re I o 1 1 o n a h i pa 

(e.g., spotiol orrongement, verticol extent of 
the observation cell) and the attributes os 
models should then serve to enhonce the 
sensed observotions for collecting informotion obout 
exomple, this opprooch should leod to on improvement 


between londscope factors 
londscope elements within 
measured in the doto. These 
utility of remotely 
the surface. For 


i ty to use the doto to identify ond map nominol closses by 
0 physical basis for explaining ond predicting the degree to 


in our obi 
providing 

which ony ground cotegory Is distinguishable 
londscope. This knowledge should also serve 
of remotely sensed doto oppropriote to 
odditlon, if the doto provide consistent 
ottributes, nomothetic solutions to the 
should be possible which will 
geographic location, with a m 


from others present in the 
os 0 guide in the selection 
0 particular 


problem. In 
meosurements of londscope 
inttirpretotion of the doto 
ollow onolysls t'f the doto from ony 
nimum of ancillary informotion. 


The potential of remote sensing os a new source of geogrophic 
knowledge obout the environment hos not been fully exploited (Bowden, 
1877). We expect thot the ottributes measurable with remotely sensed 
doto will aerve os new indicotors of londscope conditions. For exompte, 
if these obaervotions provide o meosure of the degree to which o 
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londBcope 7» phot9aynthotlco It y octtvo. w« Hove o new ffloons to study the 
biophyaicol functioning of landscapes In relotlon to noturol and 
cutturol foctors. This Intogroted expression of londssope dynomlcs 
should be highly omenoble to modeling, ond model refinement should be 
focllltoted oy the ovolloblllty of o large doto base ogolnst which to 
test model predictions, The results of this reseorch could leod to 
significant new geographic ioneepts thut would osslst In environmental 
analysis ond resources ossessment as well os contribute to ndvonced 
understanding of londscopeti on a globol basis, 

IMPLICATIONS FOR REMOTE SENSING EDUCATION 

Within mony geogrophy doportmonts remote sensing Is viewed os o 
mere "technique" a student should leorn In order to corry out "true" 
gcogrophio reseorch. This view inhibits both students ond faculty from 
I n ve 3 t I g 0 t i on of remotely sensed doto os o new source of geogrophlc 
knowledge thot moy o I te r our un de r s to n d i n g of the earth. The tendency 
hoa been for geogrophers to occept these new doto ond onolysis 
techniques from engineers and mo themo t1 c I ons without questioning the 
accompanying promises. This "block-box" opprooch hos hindered 

geographic opplicotlons of the new remotely sensed data and hos limited 
the geographer's contribution to further development rf remote sensing 
observotion systems. 

We suggest that geogrophers occept their Inherent responsibility to 
contribute to the development of remote sensing through pursuit of basic 
research along the lines we hove proposed. This reseorch can be 
encouroged, particularly among students, by demonstroting tne links 
between geogrophic theory ond remotely sensed obse r vo tl ons , encouraging 
0 hflolthy skepticism concerning our current un de r s ton d I ng of these doto, 
ond suggesting possible ovenues of research which mo)' Improve our 
understanding. The Incorporation of the fromework of Inquiry proposed 
here into current geographic remote sensing research ond education 
presents o chollenge. Rising to this challenge will, by bringing the 
full weight of the geographic perspective to beor on these new 
observations, contribute to the realization of the inherent value of 
contemporory ond future remote sensing systems. At the some time, 

pursuit of onswers to questions such os those we have posed should 
enhance our u n d a r s t o n d I n g of landscapes. 
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